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ABSTRACT 

A liquid crystal thermographic technique has been devel- 
oped which provides an excellent means of obtaining both 
qualitative and quantitative heat transfer information on 
heated objects placed in forced convection environments. 

Circumferential variation of the Nusselt number on a uni- 
formly heated right circular cylinder cooled by forced con- 
vection was obtained for Reynolds numbers varying from 38,000 
to 148,000. The results compare within the experimental un- 
certainty in forward stagnation regions with the theory of 
Schuh. Beyond approximately 30°, the results diverge from 
theory but are consistent with the work of other investiga- 
tors . 
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I . INTRODUCTION 



The purpose of this investigation was to develop a tech- 
nique for the visual determination of both qualitative and 
quantitative heat transfer and fluid flow information on 
heated objects placed in forced convection environments. 
Cholesteric liquid crystals, a material that exhibits bril- 
liant changes in color over discrete, reproducible tempera- 
ture bands, were used as the temperature sensor in the 
technique. The liquid crystals allow one to visually observe 
select isotherms and, further, can be used to infer the loca- 
tion of points of flow separation and boundary layer reattach- 
ment. The liquid crystals also give a dramatic indication 
of the influence of turbulence on surface temperature. 

The technique has been used to determine the circumferen- 
tial variation of the Nusselt number on a uniformly heated 
right circular cylinder cooled by forced convection in a 
cross flow of air. Reynolds numbers were varied from 38,000 
to 148,000 in the present investigation. This allowed the 
study of both critical and subcritical flows. The free 
stream turbulence intensity was approximately 0.5 - 0.7%. 

Data were obtained on a 4 inch diameter right circular 
cylinder constructed from .039 inch thick carbon impregnated 
paper that had a resistivity of approximately 1 ohm-in. The 
surface of the cylinder was electrically heated by passing a 
known current longitudinally through the paper. The inner 
hollow space formed by the cylinder was firmly packed with 
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glass wool to prevent heat losses into the cylinder. The 
glass wool also reinforced the cylinder and aided in resisting 
deformation due to the outer flow. 

The results obtained in the present investigation com- 
pared within the estimated experimental uncertainty in the 
forward stagnation region with a theory proposed by Schuh [1] . 
Beyond approximately 30®, the experimental values rapidly 
diverge from the theory. This is consistent with results 
found by other investigators [2,3] and is most probably ex- 
plained by the fact that the pressure distribution used in 
Schuh' s theory was taken to be the ideal, frictionless dis- 
tribution. 



II . BACKGROUND 



The use of liquid crystals in wind tunnel experiments was 
first investigated by Klein [4] in 1968. Although they showed 
great promise for obtaining heat transfer rates, several 
problems existed. The degree of accuracy was limited by the 
influences of mechanical shear, ultraviolet light, and chem- 
ical contamination. 

Some time after the completion of Klein's study, the 
National Cash Register Company (NCR) developed a process for 
encapsulating liquid crystals. McElderry [5] used these 
crystals in a study of boundary layer transition at super- 
sonic speeds, and found the encapsulated crystals to be 
relatively insensitive to shear and contamination. 
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Wirzburger (6] , in a thesis study, used liquid crystals 
to study the heat transfer characteristics of a container 
stored rocket motor placed in a hostile environment. The 
temperature distributions produced by resistively heated, 
cryogenic, and radio-frequency surgical probes were studied 
by Groff, Petrovic, and Katz [7,8,9,10,11, 12] using liquid 
crystals. 

One of the original objectives of a thesis investigation 
by Meyer [3] was to use encapsulated liquid crystals to study 
the surface temperature distribution on a uniformly heated 
cylinder cooled by a cross flow of air. This objective was 
never accomplished. The temperature distributions in his 
study were all found using thermocouples. The present inves- 
tigation is a continuation of the work initiated by Meyer. 

III. EXPERIMENTAL APPARATUS 
A. LIQUID CRYSTALS 

There are three types of liquid crystals; smectic, nem- 
atic, and cholesteric [13] . Of these three, only the choles- 
teric type was used in the present investigation. 

Microencapsulated liquid crystals provide a visual dis- 
play of temperature when applied to a black surface. The 
liquid crystal surface color will pass through the visible 
color spectrum in sequential order as a function of surface 
temperature. The temperature at which a given color appears 
is a function of the cholesterol ester formulation. By 
selecting the proper formulation the entire spectrum (red to 
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violet) can be transversed in fraction or multiple degrees 
with an accuracy of approximately 0.1*C [14]. Fergason [13, 
15-17] , in a series of articles, presents excellent discus- 
sions on the chemistry, varieties, properties, uses, and 
limitations of liquid crystals. 

The process of microencapsulation consists of enclosing 
the liquid crystals in 20-30 micron polyvinyl alcohol capsules. 
These capsules are suspended in a water based slurry. This 
process makes the liquid crystals relatively insensitive to 
mechanical shear and chemical contamination. It also reduces 
the sensitivity to viewing angle when compared to raw liquid 
crystals [5] . 

Fourteen different liquid crystal temperature formula- 
tions were used in the present study. A temperature range 
of 30*C to 49®C was covered. The liquid crystals used and 
their calibration points are listed in Table I. 

The liquid crystals were calibrated using a Rosemount 
constant temperature bath following the general procedure 
outlined by Petrovic [8]. 

All crystals were calibrated on a piece of the material 
to which they would be applied for data collection. With 
the exception of liquid crystal R-45 all liquid crystals in 
Table I were calibrated on a piece of the carbon impregnated 
paper. A piece of tape was the base material for the R-45 
calibration. 
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Table I 



Liquid Crystal Calibration Results 



Liquid Crystal Color Transition 



(NCR Designation) 


Red(®F) 


Green (®F) 


Blue(®F) 


S-30 


85.9 


86.4 


87.5 


S-32 


90.5 


91.6 


92.6 


R-33 


92.8 


93.7 


95.6 


S-34 


94.4 


95.2 


96.3 


S-36 


97.5 


98.2 


. 99.4 


R-37 


98.9 


99.7 


102.4 


S-38 


101.4 


102.4 


103.4 


S-40 


105.1 


106.1 


107.1 


R-41 


106.1 


107.4 


109.8 


S-43 


109.8 


111.0 


112.3 


S-45 


113.2 


114.0 


115.1 


R-45^ 


108.9 


109.9 


121.1 


R-49 


121.7 


122.0 


124.2 
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B, CYLINDER DESIGNS AND CONSTRUCTION 



The basic experimental procedure was first developed 
using the cylinder constructed for Meyer's work [3]. An im- 
proved procedure was established using a cylinder constructed 
from carbon impregnated paper. Final data were collected on 
a redesigned paper cylinder. 

1. Meyer ' s Cylinder 

The following is a brief description of Meyer's 
cylinder. A detailed description is contained in reference 
3 . The cylinder is shown in Figure 1 . It was constructed 
of acrylic tubing with an outer diameter of 4.47 inches. 

The cylinder was 32 inches long and spanned the wind tunnel 
test section from roof to floor. 

A surface heat flux was generated by energizing the 
0.375 inch wide, 0.003 inch thick, Nichrome ribbon which was 
helically wrapped over the middle 15 1/8 inches. Guard 
heater circuits prevented end losses and a Silastic, RTV 
foam, limited internal convective losses. 

Temperature information was obtained using thermo- 
couples welded to the inner surface of the Nichrome ribbon. 
Two pressure taps were provided for obtaining surface pres- 
sure distributions. The cylinder was mounted on a turntable 
that could be rotated from outside the wind tunnel. This 
allowed temperature and pressure information to be obtained 
at all angular locations. 

In the present set of experiments the cylinder was 
installed in the low speed Aerolab wind tunnel, located in 
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Figure 1 Meyer's Test Cylinder 



the Aeronautics Laboratory in Halligan Hall, U.S. Naval 
Postgraduate School. A schematic diagram of the wind tunnel 
is provided in Figure 2. 

The turbulence level of the wind tunnel was measured 
« 

by Meyer using a hot wire anemometer. The turbulence level 
ranged from approximately 0.5% to 0.7% for the range of 
Reynolds nvunbers studied in this experiment. 

The wind tunnel is powered by a lOOhp electric motor, 
and has a four speed transmission. The tunnel is capable of 
achieving a maximum speed of 200 mph with a clear test sec- 
tion. 

In order to obtain temperature information on Meyer ' s 
cylinder using liquid crystal, the crystals were first applied 
to tapes and the tapes were then applied to the cylinder 
surface as shown in Figure 3. 

The tapes were made by spraying Testers flat black 
paint onto Scotch double coated tape number 404. A diluted 
liquid crystal solution was then applied with a paint brush 
in two coats. Finally, the entire assembly was sealed with 
Polyurethane. The process produced tapes with an approximate 
thickness of 0.009 in. and a crystal layer of approximately 
.004 inches. 

As noted by Meyer, the surface of this cylinder had 
numerous surface irregularities caused by the difference in 
coefficients of thermal expansion between the acrylic tubing 
and the Nichrome ribbon. These irregularities can be seen 
in Figure 3, and were the primary reason for designing a 
new cylinder for continuing the investigation. 
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Figure 3. Test Cylinder with Surface Irregularities 
and the Liquid Crystal Tapes. 
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2. Preliminary Carbon Impregnated Paper Cylinder Design 

A cylinder that did not exhibit large surface irreg- 
ularities such as possessed by Meyer's cylinder was construct- 
ed from a commercially produced resistive paper known as 
Armstrong Temsheet. Temsheet is a thin, highly flexible, 
electrically resistive material. It is a carbon impregnated 
paper containing no wires or ribbons. The nominal thickness 
is 0.039 inches and the electrical resistance is approximately 
25 ohms per square. The heat that is generated when a con- 
stant current is passed through the paper is uniform to with- 
in two percent from point to point over large areas. 

A hollow cylinder with an outer diameter of 3.98 
inches and a length of 15 inches was formed from a section 
of Temsheet. 

To hold this cylinder in place two 10 inch long, 4.5 
inch outer diameter acrylic tubes were attached to the floor 
and ceiling of the 32 inch high wind tunnel test section. 

When placed inside these tubes, the Temsheet cylinder rested 
on a lip inside the lower cylinder. Two 3.9 inch diameter 
plugs threaded to a 16 inch center rod" were used to force 
the Temsheet snug against the tubing walls. The longitudi- 
nal seam of the Temsheet cylinder was sealed with tape and 
was placed at 135° from forward stagnation on the side not 
viewed through the wind tunnel window. The details of con- 
struction can be seen in Figure 4 . 

When assembled, a 12 inch long section of the Tem- 
sheet cylinder was located in the center of the wind tunnel 
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Figure 4. 
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Tube 



Cylinder with Acrylic 
Holders . 
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test section. A heat flux was generated by using a Lamba 
Regulated power supply model LK345A to pass current through 
the Temsheet. 

To insure uniform application of voltage, aluminum 
tape was used as electrodes on the inner surface of the Tem- 
sheet. The tape was placed circvimferentially on the surface 
2 inches from the top and bottom of the 12 inch test section. 
Intimate electrical contact between the tape and the Temsheet 
was established by applying a silver based conducting paint 
on the edges of the tape. This procedure resulted in a uni- 
formly heated 8 inch test section. 

To guard against internal free convection losses the 
cylinder was filled with glass wool. 

Since Temsheet is black, the liquid crystals were 
applied directly to the test surface. To insure that edge 
effects would not interfere with readings, no liquid crystals 
were placed within one inch of the electrodes. See Appendix 
B for discussion of end losses. The following liquid crystals 
bands, 3/8 inch wide and separated by 1/8 inch, were painted 
circumferentially on the surface of the Temsheet cylinder. 
Reading in order from top to bottom: R-49, S-45, S-43, R-41, 
S-40, S-38, R-37, S-36, S-34, R-33, S-32, and S-30. 

The surface of the cylinder was marked every five 
degrees to allow one to accurately locate the various iso- 
therms displayed by the liquid crystals. 



20 



3. Final Temsheet Cylinder Design 

The protruding lips on the acrylic tubes which held 
the preliminary Temsheet cylinder produced an undesirable 
secondary flow in the direction of the cylinder axis. This 
axial velocity component caused the overall flow to be highly 
three dimensional. Since only two dimensional flow was de- 
sired, a new set of supports were designed that exhibited a 
constant outer diameter to the flow. The final cylinder is 
shown in Figures 5 and 6. 

The upper and lower supports for the Temsheet test 
section were 3.96 inch diameter wooden cylinders permanently 
attached to the floor and ceiling of the wind tunnel. The 
lower cylinder had a one inch diameter hole drilled through 
its axis to allow the passage of electrical leads. 

The wooden cylinders were each 12 inches long with 
the last 2 inches on the ends turned down approximately . 04 
inches. This allowed the Temsheet to be attached to the 
cylinders with double backed tape. The outer diameter of 
the assembly was then constant over the entire wind tunnel 
height. 

The electrodes were attached as before except that 
the test section was only 6 inches long. 

The following liquid crystals were painted directly 
on the Temsheet surface; recording from top to bottom: R-49, 
S-45, S-43, S-40, S-38, S-36, S-34, and S-32. 

A second cylinder, utilizing these same bases, was 
constructed with one crystal, S-43, covering the entire sur- 
face from electrode to electrode. 
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Figure 5. Wooden Bases for Final 
Temsheet Cylinder. 
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Figure 6. Final Temsheet Cylinder 
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Finally, both cylinders were marked every five degrees 
and filled with glass wool as before. 

IV. EXPERIMENTAL PROCEDURE 

This portion of the investigation consisted of three dis- 
tinct phases. 

Initially, the cylinder used in Meyer's investigation [3] 
was used to check the liquid crystal readings against those 
of the installed thermocouples. The procedure outlined by 
Meyer was used to bring the cylinder to steady state. The 
liquid crystal temperature display was recorded and thermo- 
couple readings taken at the forward stagnation point and 
then every 15 degrees of arc. 

Some difference was expected between the two temperature 
distributions obtained. This difference was due to the loca- 
tions of the two sensors. The thermocouples were located on 
the inside of the Nichrome ribbon while the liquid crystals 
read temperature on the surface of the cylinder. Appendix C 
contains an analysis of the drop in temperature between the 
thermocouple and liquid crystal surface. 

Several items were noted during this phase of experimen- 
tation. First, the degree of agreement between thermocouple 
and liquid crystals depended on noting the point at which a 
given crystal formulation began its color transition. This 
was not suprising since calibration procedures mark only the 
beginning of color transition and a given color may be visible 
over fractions or multiple degrees Fahrenheit. 
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Secondly, the surface irregularities caused the creation 
of local hot and cold spots. This could easily be seen while 
the cylinder was wanning up or when the cylinder was energized 
with no air flow. Further, the flow pattern in the area of 
the irregularities was altered. This was visible as small 
regions of pulsating liquid crystal colors directly behind 
the irregularity. 

Additionally, the number of surface irregularities in- 
creased as the number of heating and cooling cycles on the 
cylinder were increased. 

It required approximately 12 hours for Meyer to collect 
a complete set of thermocouple data. Using the liquid crys- 
tals this time could be reduced to approximately two to four 
hours. Even this was considered too long for practical in- 
vestigations . 

It was decided to construct a new cylinder which would 
have a smooth surface and allow data to be collected within 
a reasonable time period. 

Before proceeding with a description of the experimental 
procedure using the Temsheet cylinder, a brief description 
of the temperature distribution that exists on a uniformly 
heated cylinder placed in a cross flow of air will be given. 
This description together with the accompanying Figures 7 
and 8 should prove helpful in interpreting the liquid crystal 
results. 

The flow of a real fluid is best analyzed by the changes 
occurring within the boundary layer that forms near the 
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Figure 7. Schematic of a Cylinder in Subcritical 
Flow with the Resulting Temperature and 
Heat Transfer Distributions. 
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Figure 8. Schematic Diagram of a Cylinder in 
Critical Flow with the Resulting 
Temperature and Heat Transfer Distribution 
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surface. These boundary layers are classified as either 
laminar or turbulent. The flow within a laminar layer is 
characterized by fluid particles moving in a smooth, orderly 
fashion and in layers. A turbulent boundary layer is charac- 
terized by individual fluid packets flowing in irregular 
patterns within the boundary layer. 

The heat transfer coefficients and temperature distribu- 
tions are functions of both the thickness of the boundary 
layer and the type. As the boundary layer increase in thick- 
ness with angular location, the fluid velocity is. retarded 
resulting in an increased thermal resistance to heat transfer. 
In the case of a uniformly heated cylinder the surface tem- 
perature must increase to maintain the wall heat flux at a 
constant value. The nature of turbulent flow is conductive 
to better heat transfer. 

Flow past a cylinder may be classified as either sub- 
critical or critical. 

In subcritical flow the laminar boundary layer grows 
from a minimxxm at the forward stagnation point to a maximum 
at an angular location of 80-85 degrees. At this point the 
kinetic energy of the fluid is attenuated sufficiently that 
the adverse pressure gradient present on the surface can no 
longer be overcome. The laminar boundary layer separates 
from the surface. The surface is then scrubbed by the tur- 
bulent actions of the wake. 

Figure 7 depicts a typical subcritical flow pattern and 
the resulting trends in both temperature and heat transfer 
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coefficient with angular location. Note that the maximum 
temperature exists at the separation point. This corresponds 
to the point of maximum thermal resistance, or minimum heat 
transfer coefficient. 

Critical flow is characterized by an increasing laminar 
boundary layer thickness on the forward portion of the cylin- 
der just as in the subcritical case. In the vicinity of 
80-85 degrees, however, the flow begins a transition from 
laminar to turbulent flow. The point at which transition 
begins is known as the laminar separation point. The laminar 
boundary layer separates but the kinetic energy of the fluid 
is high enough that the flow reattaches to the surface at 
some point aft of the laminar separation point. A turbulent 
boundary layer then develops and finally separates from the 
surface in the region of 110-130 degrees (the actual point 
of reattachment and ultimate separation is a function of the 
Reynolds number) . 

Figure 8 depicts the critical flow pattern and the re- 
sulting trends in both temperature and heat transfer coeffi- 
cient with angular location. 

The development of a consistent method of obtaining 
accurate heat transfer data with liquid crystals was the 
objective of the next phase of experimentation. The tempera- 
ture fields that existed on the Temsheet cylinder held by 
acrylic tubes were observed under steady state conditions in 
the wind tunnel. The best procedure for obtaining accurate 
information was developed by analyzing the visual temperature 
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fields and comparing them with those of the previously dis- 
cussed theory of angular temperature and heat transfer dis- 
tributions . 

As noted earlier, the most accurate temperatures were 
found at the beginning of color transitions on the liquid 
crystals. In areas of large temperature gradients there was 
no problem in determining the start of red, green, or blue. 
The particular liquid crystal in the temperature range at 
these points would undergo the red to green to blue transi- 
tion in a short interval. The point of blue transition was 
easily noted. 

The point of transition in the area of stagnation was 
much more difficult to observe. The temperature gradients 
in this region were very shallow. It was not uncommon for 
one liquid crystal to exhibit approximately the same color 
over 15 to 20 degrees of arc. The actual transition to this 
color was usually not apparent. The error in reading a tem- 
perature in this region could be as large as 1®F. To obtain 
a precise temperature the beginning of a transition was 
forced by adjusting the power supply until a previously un- 
changed liquid crystal just began its transition. This was 
verified by decreasing the voltage a few tenths of a volt 
and observing the liquid crystal return to its previous 
color. 

If a clear cut transition from red to green, or green to 
blue occurred in any region of shallow temperature gradients, 
it was considered precise and used. This type of transition 
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usually occurred in angular regions 15 to 60 degrees aft of 
the forward stagnation point. 

The hottest point on the cylinder surface was at separa- 
tion (Figure 7) . This point was found by adjusting voltage 
until the first point of red color was observed on an un- 
changed crystal or decreasing power on one that was already 
through its transition until the very last red point remained. 
The location was checked by repeating the procedure on a 
different liquid crystal band. The two points invariably 
fell on a vertical line. 

By combining the procedures repeatable results were ob- 
tained for a given Reynolds number. 

While experimenting with the above procedure the Reynolds 
number was maintained at approximately 50,000, clearly in the 
subcritical flow regime. However, a local cold spot not at- 
tributable to normal subcritical flow patterns existed in the 
region of 90 degrees. Small strings placed on the cylinder 
confirmed that undesirable secondary flow patterns were gen- 
erated by the lips on the acrylic tubes. Due to this flow 
anomaly a third phase of experimentation employing an im- 
proved wooden base was initiated. The objective of the third 
phase of the investigation was to gather data for comparison 
with conventional works. This was done in three steps. The 
new wooden base Temsheet cylinder was tested for two dimen- 
sional flow characteristics and edge losses. The roughness 
of the Temsheet and liquid crystal coatings were examined. 

Data for Reynolds numbers ranging from 38,000 to 148,000 were 
obtained for analysis. 
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For data collection the Temsheet cylinder with eight 
liquid crystal bands was installed in the wind tvinnel. 

The airstream thermocouple and its reference ice bath 
were inspected and temperature recorded for initial calcula- 
tion of airstream properties and velocity. 

A U-tube manometer and micromanometer used for airspeed 
indication were carefully zeroed. 

The power supply was energized and voltage adjusted to 
cause all liquid crystal bands to undergo a transition. 

This served two purposes. First, if the glass wool packing 
was not installed snuggly a point of temperature discontinu- 
ity would appear. The packing could then be adjusted prior 
to establishing air flow. Secondly, the cylinder could be 
preheated to reduce the time necessary to reach steady state. 

The wind tunnel was started and the air speed was in- 
creased using the previously calculated pressure drop on the 
U-tube manometer. 

Since the U-tube manometer was not very sensitive over 
the entire range of Reynolds numbers, especially at the lower 
Reynolds numbers, the pressure drop read on the micromanometer 
was used to recalculate air speed and Reynolds number. This 
proved to be a critical part of the experiment. 

Power was increased until the upper liquid crystal, R-49, 
began its red transition at separation. Power was continu- 
ally adjusted to maintain red at separation until steady 
state was reached. This usually occurred within fifteen 
minutes and was typified by a period of approximately five 
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or more minutes of a constant red without adjusting power. 
This was checked by decreasing power a few tenths of a volt 
euid noting the red spot disappear. 

Data were collected after restoring power to its previous 
setting. The angular location of each liquid crystal color 
transition, voltage, airstream temperature, and micromanom- 
eter height were recorded. 

For critical flows the procedure for locating separation 
was followed to locate the points of laminar separation, 
reattachment, and final separation. 

The test for edge effects and two dimensional flow char- 
acteristics consisted of placing a Temsheet cylinder coated 
with liquid crystal S-43 from electrode to electrode in the 
wind tunnel. 

The wind tunnel was started and the cylinder heated to 
steady state. The isotherms were then observed for straight- 
ness and end effects. 

Although the surface of the Temsheet was smooth to the 
touch the actual surface roughness and effects of the liquid 
crystals were not known. A piece of Temsheet coated with 
liquid crystals was examined under an electron microscope. 
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V. RESULTS 



A. TEMPERATURE DISTRIBUTION AND HEAT TRANSFER 

The nximerical data and data reduction technique are con- 
tained in Appendix A. Briefly, the local heat transfer co- 
efficient was obtained by dividing the heat generated per 
unit area by the temperature difference between the surface 
and air stream temperatures. 

h = q/AT 

The value of q was determined by dividing electrical power 
generated by the surface area 

q = (vVr) 

A 

The local convective heat transfer coefficient, h^, was then 
obtained by subtracting the radiation coefficient, h^, (see 
Appendix A) from h. 

he = h-hj. 

The Nusselt number (Nu) was calculated using the expres- 
sion 

Nu = 

k 

and the Froessling number calculated from 



Fr = 



Nu 



Data are presented graphically in this section. Figures 
of temperature distribution, Nusselt numbers, and Froessling 
numbers are included. 
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Figure 9 shows a sketch of a typical subcritical flow 
temperature distribution and a representation of the isotherms 
seen on the liquid crystal bands on the cylinder surface. 

Critical flow was noted at a Reynolds number of approxi- 
mately 120,000. The shape of a typical critical flow temper- 
ature distribution was shown in Figure 8. Figure 10 is a 
color photo of the actual temperature distribution on the 
cylinder surface for this critical flow condition. Note that 
on the upper liquid crystal band, R-49, the presence of the 
separation region is characterized by a maximum temperature 
at 115° (denoted by a red patch) with lower temperatures on 
either side (denoted by black) . The blue region forward of 
this region indicate higher temperatures and bracket the 
laminar separation point. 

The Nusselt numbers for several Reynolds numbers are 
plotted in Figure 11. Several items should be noted from 
these curves. 

First, the general shape of the curves is consistent 
with the previous work of Giedt and Meyer [2,3] in both crit- 
ical and subcritical flows. 

The points of minimum heat transfer coefficients occur 
in the region of 80 to 90®. However, unlike the work of 
Giedt, these points move forward on the cylinder with in- 
creasing Reynolds numbers. This may be due to a slight 
deformation of the cylinder as velocity increases. Deforma- 
tion was not visible to the naked eye at Reynolds numbers 
below 150,000. It was interesting to note that the isotherms 
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Cylinder Surface Liquid Crystal Bands 




Figure 9. Liquid Crystal Temperature Field and Resulting 

Temperature Distribution on the Cylinder Surface. 
Subcritical flow conditions indicated. 
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Figure 10. Liquid Crystal Temperature 
Field in Critical Flow. 
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O Re = 38,000 
Q Re = 88,000 




Figure 11. Heat Transfer Results at Reynolds Numbers 
of 38,000, 88,000, and 148,000. 
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depicted on the cylinder surface were no longer vertical 
when deformation was visible. This can be seen in Figure 12. 

The Froessling numbers for the laminar flow regions are 
plotted in Figure 13 for comparison with the theory of Schuh 
[1] . Comparison is within experimental uncertainty near 
stagnation and then drops off near separation. The latter 
trend is consistent with the work of Giedt and is most prob- 
ably explained by the fact that Schuh used an ideal pressure 
distribution as opposed to the distribution encountered in 
viscous flows. 

B . EDGE EFFECTS 

Figure 14 showed the isotherms on the cylinder coated 
with one liquid crystal, S-43. They clearly indicate that 
heat is lost near the top and bottom edges of the coated 
region. 

The edge effects are indicated by a slight curvature of 
the isotherms near the electrodes. This is expected from 
the model developed in Appendix B and was the reason for not 
taking data within an inch of the electrodes. 

C. SURFACE ROUGHNESS OF TEMSHEET 

Figures 15 and 16 show the Temsheet liquid crystal inter- 
face magnified 121 and 27 times respectively. The point to 
note is that the liquid crystals fill-in irregularities in 
the surface and appear to be 20-30 micron spheres as speci- 
fied by the manufacturer. 
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Figure 12. Cylinder Deformation at Reynolds 
Numbers Greater than 150,000. 
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X Re 
O Re 
A Re 
□ Re 
• Re 
■ Re 



38.000 

48.000 

68.000 

87.000 

108.000 

121,000 




Angular Location 

Figure 13. Comparison of Experimental Results of the 
Present Investigation with the Theoretical 
Solution of Schuh. 
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Figure 14. Isotherms Showing End Losses and 
Two Dimensional Flow Pattern. 
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Figure 15. Temsheet-Liquid Crystal Surface 

at 121X. 




Figure 16. Temsheet-Liquid Crystal Surface 

at 27X. 
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Page and Warsap in reference 18 predict a transition to 
critical flow at a Reynolds number of approximately 148,000 
for flow over a cylinder with a surface roughness equivalent 
to the roughness of the crystal coated Temsheet (k/D=2xio~‘* ) ; 
k being the critical roughness height and D the diameter of 
the cylinder. The actual transition occurred at approxi- 
mately 120,000. 



VI. CONCLUSIONS AND RECOMMENDATIONS 

The liquid crystal thermographic technique developed in 
this investigation provides an excellent means of obtaining 
both qualitative and quantitative heat transfer information 
on heated objects placed in forced convection environments. 
Using the technique it was possible to quickly and easily 
obtain information on the variation of the Nusselt number 
around the circumference of a uniformly heated right circular 
cylinder placed in a cross flow of air. The technique also 
allowed one to visually observe the effects of flow separa- 
tion, the turbulent boundary layer, and the turbulent wake 
on the surface temperature of the cylinder. Colored movies 
taken of a cylinder coated with a single liquid crystal are 
especially vivid in their display of the influence of the 
turbulent wake on the cylinder surface temperature. In the 
wake region the crystals alternately dim and glow in response 
to the "scrubbing" action caused by random bursts of cool 
fluid impacting on the surface. 
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The ability to visually observe turbulent flow pattern 
effects on surface temperatures suggests an excellent means 
for studying the influence of free stream turbulence on the 
heat transfer rates of heated objects. However, in order 
that such a study be conducted in a quantitative manner, the 
thermal response time of the liquid crystals must first be 
determined. Parker [19], using a capacitor discharge tech- 
nique, found that a 0.001 inch film of liquid crystals 
coated on a thin stainless steel foil responded in approxi- 
mately 0.036 seconds to a step change in foil temperature. 

It is possible that a thin coat of liquid crystals placed on 
a material such as Temsheet would respond even faster than 
this due to the fact that the crystals are partially ab- 
sorbed by the Temsheet. The ideal situation would be if the 
liquid crystals responded at the same rate as the Temsheet 
itself. An experiment, similar to the one conducted by Parker, 
needs to be conducted to determine response time information 
on the liquid crystals. 

As a final recommendation, the phenomeron shown in Figure 
17 is offered as a possible topic for future investigation. 
During the final phase of collecting wind tunnel data for 
the present investigation, it was noted that the cylinder 
coated with liquid crystal S-43 displayed alternate hot and 
cold spots along the separation line. These spots were uni- 
formly spaced and seemed to be caused by some flow phenomenon, 
perhaps a series of vortices such as one observes in flow in 
a curved channel. Whatever the cause, the hot and cold spots 
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Figure 17. Hot and Cold Spots on the Cylinder 
Surface along Separation Line. 
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definitely existed as is readly apparent in Figure 17. Pre- 
cise measurements were not taken at the time the phenomenon 
was observed (the Reynolds number was approximately 75,000) 
and no explanation is offered for its existence. It should 
be noted, however, that such a phenomenon may well have gone 
undetected if thermocouples were used as the temperature 
sensor. Additional research, perhaps using liquid crystals 
for temperature sensing and smoke for flow visualization, 
needs to be conducted to explain the observed phenomenon. 
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APPENDIX A 



DATA AND DATA REDUCTION 

In order to experimentally determine the Nusselt number 
as a function of angular location on the right circular 
cylinder under investigation, a relationship between the con- 
vective heat transfer coefficient, h^, and the surface heat 
flux produced by the Joulean heating effect in the Temsheet, 
V^/RA, was needed. This relationship is easily developed by 
performing a simple energy balance on an elemental volume, 
see Giedt [2] or Meyer [3] . The result is 

he = RA r^^ de^ - h^ 

(T-T„) 

where 

h^ = Convective heat transfer coefficient 
V = Voltage impressed across the test section 
R = Electrical resistance of test section 
A = Area of test section 
t = Thickness of test section 

k = Thermal conductivity of test section material 

r^ = Radius of the cylinder 

T = Temperature at angular location 0 

0 = Angular location 

Tjo = Air stream temperature 

h_ = Radiation heat transfer coefficient = oF. ,{T+T }{T^+T 
Fi _2 = Radiation exchange factor between cylinder and sur- 
roundings 
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a = Stefan-Boltzman Constant 

An exact value for thermal conductivity of the Temsheet 
material was not available, but assuming that its value is 
approximately the same as other carbon paper products, the 
product kt is calculated to be .000195 BTU/hr-®F. Since this 
value is so low it was felt that the conduction term in the 
expression for convective coefficient could be neglected. 

The expression used for the convective coefficient then takes 
the final form; 



he = 






- h. 



RA (T-T«) 



The Nusselt number is: 

Nu = ^h^rQ 
^air 

The Froessling number is: 

Fr = Nu/ 11^ 

where 

Re = Reynolds number = V<joD/v 
V<» = Air velocity 
V = Air kinematic viscosity 

A sample calculation is provided for illustration using 
the following values from the data at a Reynolds mamber 
47,818: 



0 = 60* 
e = .9 

a = .1714 X 10“®BTU/hr-ft^-®R“ 
To = 62.7°F 
T = 96.3°F 
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Vco = 22.5 ft/sec 

V =23.7 Volts 

R = 12.0 n 

D/K = 22 hr-®F-ft2 
BTU 

Re = 47,818 

h = (23.7V) ^ 

(12.0A) (.52ft^) (33. 6 “F) ( . 293WATTS/BTU/hr 

h =9.2 BTU 

hr-ft^'-^F 

h = .9(.1714 X 10-®BTU/hr-ft^-*’R'*) (556.3+522.7) 

(556. 3^+552. 7^) 

hy. = .97 BTU 

hr-ft2-*F 

h =8.2 BTU 

hr-ft2-®F 

Nu = (8.2 BTU ) (22 hr-*F-ft^) 

hr-ft2-*F BTU 

Nu = 180.4 

Fr = 180.4 = .82 

“^47,818 

The data sheets for all seven Reynolds numbers investi- 
gated follow. 
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DATA 



Air Temp. = ||c 

V ft 

Air Speed = •/? onHaO 
* ft/sec 
Re = 36500 

Correction Factor = /■0^(>y 

Recorr* 3 
Voltage = 3^6" 4 V 
Heat Flux BTU/hr.ft2 
Film Temp = ^6) *F; 
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DATA 



Air Temp. -L>2.1 ®F; = 2010 

Air Speed = d23 anH20 
ft/sec 

Re = A(s>,S7S’ 

Correction Factor - /,OZl>2 

Re = 
corr ^ 

Voltage = y 

Heat Flux -20Q ^T^^/hr.ft^ 

Film Temp = 85" T; 



sec 

ft 



Date: 

Time: /OOO 

Crystal reader: FIELD 
Data sheet: CooPt/i. 

Cylinder Dia. = O. SZo /-r 
Cylinder Area = O. 
Resistance = IZ.Qjl 
Crystals: 
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DATA 



Air Temp. = *F; -?.o7q 

Air Speed - . C, crniHaO 
» 3?. ft/sec 
Re = Cc>j^^O 

Correction Factor = 



Re 



corr 



^5. oo$ 



^ Date: J-7.(o-7/ 

Time: 

Crystal reader: 

Data sheet: Coop^^ 
Cylinder Dia. = . 38 
Cylinder Area = ° 
Resistance = /s.o-^ 



Voltage = V ' uf sjy 

Heat Flux • 

S3y 

Film Temp = 
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DATA 



7 3 

V 



Air Temp, = Q't-l *F; 

Air Speed = cmH20 

® ft/sec 

Rg = /dS‘,0<a0 

Correction Factor = /<o^l 

Re,_ = /o? sc^S' 

corr ' 

Voltage = V 

Heat Flux BTU/hr.ft2 

Film Temp = ®F; ^ 



sec 

ft 



Date: 

Time: /3 as“ 

Crystal reader: 

Data sheet: Coo/>&^ 
Cylinder Dia. = .3iPt 
Cylinder Area = 
Resistance = /^-c> 
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/749 

5<^5" 

^'-lO 



S 3y 
63 c. 
55/ 
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APPENDIX B 



EDGE EFFECTS 

In the set of experiments involving the cylinders con- 
structed from Temsheet, the top and bottom edges of the test 
section were not guarded against heat losses. An analysis 
was performed and it was found that edge effects did not in- 
fluence the temperature field beyond a distance of approxi- 
mately 1/4 inch. The details of the analysis are presented 
in this appendix. 




The axial variation in the temperature field due to end 
losses was estimated by treating the Temsheet as a thin fin. 
Symmetry about the centerline was assumed. Temperature 
variations in the circumferential direction were neglected. 

A "worst case" situation was studied by assuming the temper- 
ature at y=L equaled the free stream air temperature, T^. 
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Due to syiranetry, the centerline represented an adiabatic sur- 
face. The governing equation for this situation is: 

+ i" ' - h , 0 (1) 

dy * k kt ~ 



where 

T = Temperature at location y 

' = Volumetric heat generation rate = V^/2RLWt 

V = Voltage impressed across the electrodes 

R = Electrical resistance of material between electrodes 

h = h^+h_ = Surface heat transfer coefficient 
c r 

k = Thermal conductivity of test section material 
t = Thickness of test section material 
Too = Air stream temperature 

The boundary conditions on the problem are: 



at y = 0, 


dT/dy = 0 


(2) 


at y = L, 


T = T 
•** ***00 


(3) 



It should also be noted that the surface of the Temsheet 

that was in contact with the glass wool was assumed to be 

perfectly insulated. The solution to equation (1) that 

satisfies boundary conditions (2) and (3) is 

T-Too = _ Cosh (my)) 

T^eo Cosh(mL) (4) 

where 

• f t f 

(Tq-Too) = Centerline temperature excess = q t 




Equation (4) was used to determine the minimum distance 
from the electrodes that data could be taken and not be 
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affected to any substantial degree by end losses. As a guide- 
line, equation (4) was used to determine the y location for 
which T-T /Tq-Too = 0.99. Using representative values of; 

2L — 6 inches 
h = 5 BTU/hr-ft^-®F 
k = 0.06 BTU/hr-ft-^F 
t = 0.00325 ft 

This distance is found to be: 
y = 2.78 inches. 

It was concluded that if the electrodes were spaced 6 
inches apart, using only the center 4 inches would yield 
temperature information that was essentially uninfluenced 
by edge effects. Figure 14 in the body of the Thesis sup- 
ports this claim. 
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APPENDIX C 



COMPARISON OF THERMOCOUPLES AND LIQUID CRYSTALS 

In the set of experiments that employed the Nichrome 
wrapped cylinder, a comparison was made between the tempera- 
ture determined with the thermocouples and the temperature 
as determined visually using the liquid crystal tapes. This 
comparison was made to serve as a check on the accuracy with 
which one can determine temperature using liquid crystals. 

In comparing the two temperature sensing techniques it must 
be noted that the thermocouples were located on the inside 
of the Nichrome ribbon whereas the liquid crystals were lo- 
cated on the outside of the tape. As such, a thermal resis- 
tance exists between the two locations. This will result 
in a temperature drop. An elementary heat transfer analysis 
on the composite system as shown sketched below yields the 
following expression for this temperature drop 



1 


n 


kj 




Lt 1 2 

1 


Tape Composite 




r V 


Nichrome i 



y//v//////iiiii ininiininni 






The subscripts "1" refers to the Nichrome and the sub- 
script "2” refers to the tape. The volumetric heat 
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generation term, 
the ribbon as q ' ' ' = 
where 



is related to the voltage drop across 



V^/RAt 



R = Electrical resistance of the ribbon 
A = Total surface area of the ribbon 
t = Thiclcness of the ribbon 
V = Voltage 

The following were taken as typical thicknesses and properties 
tj = .00025 ft 
tz = .00075 ft 
kj = 5.5 BTU/hr-ft-*F 
kj = .12 BTU/hr-ft-®F 

Using these values the expression for AT becomes 
AT = 2.07 X 10"®q''' ["F] 

where q' ' ' is expressed in BTU/hr-ft^ 

The above expression was then used to evaluate the re- 
sults of the test. Several items where noted. The thermo- 
couple readings should have been higher than those of the 
liquid crystals. This was not always the case. Careful 
analysis of the surface showed that the thermocouples were 
only under tape R-45. The liquid crystal temperature dis- 
tribution was taken from several tapes. This meant that the 
temperatures compared were at the same angular location but 
not necessarily the same vertical location on the cylinder 
surface. 

It should be noted that the location of the thermocouples 
under tape R-45 was easily found. The thermocouples acted 
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as a heat sink and could easily be found during warm-up of 
the cylinder. One showed up as a pin point size cool spot 
at an angular location of 90®. 

Examination of the surface of the cylinder showed that 
the thermocouple was located on a smooth surface whereas the 
surface under the tapes directly above or below it could 
have irregular shape because of the Nichrome ribbon irregu- 
larities. The surface irregularities produced hot or cold 
spots and accounted for the departure from expected results. 

Since thermocouple position on tape R-45 was known and 
the surface there was smooth, it was decided to compare tem- 
perature readings at that point. The readings compared 
within the calculated differences. Table II shows this 
comparison. 

Table II 

Comparison of Thermocouple and 
Liquid Crystal Temperature Readings 

Temperatures (®F) Temperature Differences (®F) 



Liquid Crystals 


Thermocouple 


Maximum Predicted 


Actual 


108.9 (Red) 


109.8 


2.7 


.9 


109.9 (Green) 


112.5 


2.9 


2.6 


112.1 (Blue) 


113.5 


3.2 


1.4 
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APPENDIX D 



UNCERTAINTY ANALYSIS 

The degree of uncertainty for the final results was cal- 
culated using the method of Kline and McClintock described 
in Ref. 20. 

The measured variables which were the origins of uncer- 
tainty were voltage, resistance, surface temperature, angular 
location, micromanometer pressure drop, cylinder physical 
measurements, and the properties of air. 

The error in determining the properties of air was 
negligible. 



THE HEAT TRANSFER COEFFICIENT 



The heat transfer coefficient was calculated from 



h = vVR 

AAT 




For the following experimental run at a Reynolds number of 



47,818 



V = 


27.6 ± 


.3 


VOLTS 


R = 


12.0 ± 


.3 


OHMS 


AT = 


39.7 ± 


.9 


OF 


A = 


.52 ± . 


01 


ft 



(20 to 1) 
(20 to 1) 
(20 to 1) 
(20 to 1) 



64 





I 



• . ^ 



I 




Substituting 







/ 2 ^ • 3 . 
^ 2776 ^ 



(— ) 



( 



.9 



39.7 



) 



.01 

52 



2 



= .045 
”h 



This represents the uncertainty at the forward stagnation 
without regard to the uncertainty in determining the angular 
location. It is estimated that the uncertainty in reading 
angular location in this region would be ±5* (20 to 1) . 

The uncertainty in the heat transfer coefficient, at 
separation was approximately 4.5 percent with an angular 
uncertainty of ±1* (20 to 1) . 

The uncertainty in the wake was on the order of 8 percent 
with an angular uncertainty of ±5® (20 to 1) . 



THE NUSSELT NUMBER 

The Nusselt number was calculated from the equation 

Nu = h£ 

K 

The uncertainty equation would be 



“Nu + (“Dv ^ ^ 

Nu y h^ ^ 

Since the uncertainty in D and K was considered negligible 



%u = “h 
Nu h 

or the same as for heat transfer coefficients. 
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THE REYNOLDS NUMBER 



Reynolds number was calculated from the equation 

Re = VD 
V 

and its uncertainty is the same as free stream velocity, V. 



V 




2 Ap 
P 



“Re = “V =*lf ,.5“Apv ^ = .5 ,“Ap, 

-rS ^ (I ‘ 

Since Ap is directly related to manometer reading Ah 



“Re = .5 “Ah 
Re h 

For a Reynolds number of 37,885 



“Re = .5 ..01 = .025 

Re \19^ 

For a Reynolds number of 148,000 



“Re = .5 ..01. = .002 
Re 4.8^ 

THE FROESSLING NUMBER 



The Froessling Number was calculated from Fr = ^ 

iPe 



Then, 

“Fr =~\j “nu ^ + ,“Re ^ 
Fr U' Nu' ^Re' 
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Substituting 



^Fr 

Fr 



”^(.045) ^ + (.017) 2 = .05 



This is for the Reynolds number 47,818 used previously. It 
would decrease with increasing Reynolds nvimbers. 
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